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ABSTRACT 


The  basic  differential  equations  for  the  electron  density  in 
an  electron  beam  stabilized  discharge  are  formulated  and 
a  discussion  of  the  relevant  physical  processes  is  presented. 
Rapid  charge  transfer  among  the  positive  ions  is  shown  to 
result  in  a  considerable  simplification  of  the  basic  equations. 
An  analytic  solution  is  given  which,  in  an  approximate 
manner,  includes  the  effect  of  oxygen  attachment.  A 
comparison  between  calculated  electron  densities  and 
preliminary  measurements  is  shown  to  be  in  reasonable 
agreement. 


1.0  INTRODUCTION 


Molecular  gas  lasers,  in  which  the  primary  pumping  mechanism  is  electron 
impact  on  the  gas  molecules,  have  shown  very  high  efficiencies;  for  example 
up  to  47%  for  CO  and  33%  for  CO^.  Because  of  the  high  efficiency  and 
relative  simplicity  of  these  electrically  pumped  lasers,  considerable  effort 
is  being  directed  to  the  problem  of  scaling  up  the  outputs  of  electric  dis¬ 
charge  lasers  by  new  methods  of  maintaining  uniform  discharges  in  high 
pressure  gas  mixtures.  One  such  new  method  is  the  electron  beam  stabilized 
discharge.  ’  ’  ’  The  calculation  of  the  electron  density  of  the  plasma  in 
such  a  discharge  is  the  subject  of  this  repart. 

An  electron  beam  stabilized  plasma  differs  considerably  from  a  self- 
sustaining  discharge  in  that  the  energy  for  the  ionization  of  the  gas  molecules 
comes  from  a  Ugh  voltage  electron  beam  directed  into  the  laser  chamber 
through  a'  thin  metalic  foil  from  a  high  vacuum  region  (see  Figure  1).  In  the 
more  common  self-sustaining  discharge,  e.  g. ,  glow  or  radio  frequency 
discharges,  the  ionization  is  produced  by  electrons  generated  in  t;ie  plasma 
which  po  sess  sufficient  energy.  At  high  pressure,  for  exam  '  t  - 
atmosphere,  and  at  discharge  times  longer  than  about  1  /is,  such  plasmas 
with  metalic  electrodes  often  produce  low  resistance  local  arcs  which  in 
turn  cause  very  inhomogeneous  laser  pumping  and  unwanted  heating  of  the 
active  molecules.  Since  the  ionization  in  an  electron  beam  plasma  is 
externally  controlled  by  the  energetic  electron  stream  ,one  need  not 
"heat"  the  electrons  in  the  plasma  to  the  point  of  causing  significant 
ionization  and  thus  the  problem  of  arcing  in  the  laser  chamber  is  greatly 
reduced.  Note  that  we  distinguish  between  th<j  tenuous  stream  of  high 
energy  electrons  of  the  E-beam  and  the  low  energy  (plasma)  electrons 
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Figure  1.  Cross  section  through  the  discharge  region. 


produced  in  the  laser  chamber  by  ionizing  collisions  between  molecules  and 

the  fast  beam  particles,  (The  electron  density  of  the  beam  may  be  about 

6  3  12  3 

10  /cm  and  of  the  low  energy  electrons  about  10  /cm  ). 

Fast  laser  pumping  rates  are  obtained  if  the  average  energy  of  the  plasma 
electron  is  about  1  eV  (see  References  7  and  8  for  CO  and  CO2).  This 
average  energy  is  supplied  by  passing  a  current  between  two  su stainer 
electrodes  located  in  the  plasma  region.  One  electrode,  which  does  not 
block  the  electron  beam  is  located  near  the  foil  and  is  generally  at  ground 
potential  and  the  other  electrode  is  maintained  at  the  sustainer  potential  Vfi. 

The  laser  pumping  rate  is  proportional  to  the  number  density  of  the  plasma 
electrons,  n  ,  and  it  is  the  purpose  of  this  report  to  discuss  the  physical 
processes  which  control  the  electron  density  i  a  mcltigas  E-beam  plasma, 
and  to  calculate  n  for  some  typical  operating  conditions.  In  Section  2.0 
of  this  report  the  basic  equations  for  the  electron  density  are  formulated 
and  ng  is  calculated  for  some  exemplary  cases.  Also  a  comparison 
between  experiment  and  calculation  is  made.  A  simple  discussion  of  the 
basic  physical  processes  of  electron  generation  and  loss  is  presented  in 
Section  3.0,  The  rates  for  these  processes  are  also  calculated. 
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2.0  ELECTRON  DENSITY  CALCULATIONS 

Molecular  ga/.  laser  plasmas  generally  contain  a  mixture  of  gases,  their 

positive  ions  (and  possibly  negative  ions)  and  electrons.  In  electron  beam 

stabilized  plasmas  the  total  pressure  may  bo  about  one  atmosphere 
19  3 

(^10  particles/cm  )  or  higher.  The  number  densities  of  the  ions  or 

14  3 

electrons  is  usually  less  than  10  /cm  for  a  total  pressure  of  one 
atmosphere  and  thus  is  small  compared  to  the  neutral  density.  Since 
electron-neutral  collisions  are  esponsible  for  the  laser  pumping,  with 
ion-neutral  processes  relatively  unimportant,  the  electron  density  is 
of  prime  importance.  Although  the  concentration  of  the  various  ion  types 
is,  per  se,  not  of  interest,  we  shall  nonetheless  have  to  consider  the  various 
ion  processes  since  the  rate  of  electron  loss  depends  upon  the  particular  ion 
species  present  in  the  discharge. 

In  Section  2. 1  the  basic  set  of  equations  which  should  be  solved  to  calculate  the 
electron  density  is  given.  For  each  constituent  of  the  gas  mixture,  there 
is  one  first  order  nonlinear  differential  equation,  and  in  general,  the  result¬ 
ing  set  of  coupled  differential  equations  must  be  solved  numerically.  When 
a  chemically  pure  gas  is  considered,  one  can  write  an  analytic  solution  for  ne 
as  a  function  of  time  (Section  2.2).  One  can  also  solve  analytically  an  approxi¬ 
mate  equation  for  a  single  constituent  gas  with  a  small  amount  of  oxygen 
(Section  2.  2).  Solutions  for  a  t nu It i component  gas,  on  the  other  hand,  are 
more  complex.  A  discussion  of  the  dominant  mechanisms  in  typical  mixtures 
is  given  in  Section  2.  3.  Computer  solutions  for  various  gas  mixtures  and 
operating  conditions  are  presented  in  Section  2.  4,  and  in  Section  2.  5  some 
experimental  measurements  are  compared  with  calculated  values.  \ 

I 
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The  rate  constants  used  in  the  calculations  of  electron  density  are  identified  j 

and  reviewed  in  Section  3.  0.  Since  we  will  use  these  constants  in  the  next  j 

I 

seetion,  it  is  useful  to  briefly  describe  them  here.  | 
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Electron  Generation,  The  only  electron  generation  process  of 
importance  is  that  of  electron  neutral  collisions.  We  choose  to  divide 
the  electrons  into  two  classes.  Class  one  includes  the  electrons  of  the 
high  voltage  beam  and  the  fast  secondaries  which  the  high  voltage  beam 
produces.  Class  two  includes  the  low  energy  electrons  produced  by  the 
ionization  which  reside  in  the  plasma  at  an  energy  of  approximately  1  eV. 

Electron  Attachment.  The  attachment  processes  are  twoJbody  and  three 
body  processes.  A  two  body  process,  for  example,  is  O2  +  t.  +  O,  and  a 
three  body  process  may  be  illustrated  by  using  nitrogen:  N2  -5-  +  e  ♦  N2  +  O 

Recombination.  The  main  electron-ion  recombination  processes  are  also 
cwo  body  and  three  body  events,  illustrated  with  argon: 

Ar*  +  e  «•>  2Ar 
Ar  +e^Ar 
Ar*  +  Ar  +  e*+2Ar 

Ion-ion  recombination  is  the  main  neutralization  process  for  negative 
oxygen  ions.  The  recombination  coefficient  for  the  latter  ia  about 
10  cm  / s. 

Charge  Transfer:  The  charge  transfer  processes  are  very  important 
in  determining  th'  lifetime  of  the  various  species  of  ions.  Rate  constants 
given  in  the  works  of  Ferguson,  Schmeltekopf,  Fensenfeld  ct  al  (cited 
in  Section  3.  S)  indicate  that  there  is  rapid  charge  transfer  from  helium 
and  nitrogen  ions  to  carbon  monoxide  and  caroon  dioxide  ions. 


2. 1  Basic  Electron  Density  Equations  for  a  Multicomponent  Gas.  To 
formulate  a  basic  set  of  equations  for  a  multicomponent  gas  we  shall 
assume  qua  si -neutrality  and  homogeneity  Quasi-neutrality  implies 
near  equality  between  free  positive  and  vjgative  charges,  and  homogeneity 
implies  that  the  temporal  rate  of  char  :  of  ion  number  density  is  governed 
only  by  ion  loss  and  generation  rates.  Let  n(j^)  represent  the 
number  density  of  positive  ions  of  the  jth  species,  n  the  electron  density, 
and  n(0~  )  the  number  density  of  charged  oxygen  molecules  as  the  only 
negative  ions  present.  From  the  quasi -neutrality  requirements,  we  can 
write 

ne  =  ■?  n(j+)  -  n(0“)  (la) 

and  from  the  homogeneity  concept  we  write  (for  each  particle  species) 

dn(j+)/dt  =  G(j+)  -  L(j+)  (lb) 

and 

dn(0“)/dt  =  G(0“)  -  L(0^)  (1c) 

where  the  G  and  L  terms  refer  to  the  generation  and  loss  rates  for 
each  type  of  ion. 

There  are  many  possible  reactions  in  a  typical  laser  gas  which  must  be 
included  in  the  loss  and  generation  terms  of  Eqs.  (lb)  and  (lc).  In  order 
to  obtain  a  tractable  set  of  equations  we  shall  confine  our  attention  to 
those  reactions  listed  below  and  discussed  in  Section  3,  5. 


For  a  pure  gas,  for  example  helium,  the  dominant  reactions  for  our 
purposes  are: 


Ionization: 


e  ■§■  He  -¥■  He  *  2e 


Molecule  Formation:  He+  +  He  -*•  He* 


Recombination : 


He^  +  e-*2He. 


For  a  gas  mixture,  for  example  helium,  nitrogen  and  carbon  monoxide, 
the  dominant  reactions  are: 


Ionization: 


e  +  He  -*■  He  4-  2e, 
e  +  N*  -f  2e, 
e  +  CO  CO  +  2e. 


Charge  Exchange: 


Cluster  Formation: 


Recombination : 


He*  +  N2  (or  CO)  —N*  (or  CO*) 
N*  -?•  CO  -*-n2  +  CO+ 

CO*  +  CO~>>  (CO)*  , 

N2  +  N2^  N4' 

(CO)*  +  e  2CO, 

4- 

N.  +  e-*  2N,. 

4  c 


For  a  gas  mixture  containing  a  small  amount  of  oxygen,  in  addition  to  the 
above,  we  have: 

Oxygen  Attachment:  0,  +  e  -*0  +  0~, 

Cf 

0?  +  e  +  M  O,  +  M, 


Recombination : 


(O)  +  (CO)  ^0  +(CO>, 
x  y  x  y 


where  at  this  stage  of  the  investigation  the  specific  ions  and  molecules, 
designated  by  the  subscripts  x  and  y,  cannot  clearly  be  identified  and  M 
is  any  atom  or  molecule. 
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The  charge  exchange  processes,  schematically  indicated  above,  enter  into 
the  ion  density  rate  equations  (lb)  as  a  loss  for  the  initial  ion  species  and 
a  generation  for  the  resultant  species. 


The  basic  equations  (la,  b  and  c)  can  be  simplified  by  considering  the 

approximate  rates  of  the  various  processes.  In  particular,  the  average 

lifetime  for  helium  or  nitrogen  ions,  considering  charge  transfer  processes, 

is  in  the  order  of  10  ^sec,  and  the  lifetime,  determined  by  recombination 
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processes,  is  in  the  order  of  10  to  10  sec.  Since  we  are  not  interested 
in  the  ion  density  but  rather  electron  density,  and  since  the  free  electrons 
are  not  involved  in  the  charge  transfer  process,  we  can  simplify  the  basic 
equations.  We  sum  up  the  effect  on  the  basic  equations  of  the  dominant 
reactions  given  above  as  follows  (neglecting  for  the  moment,  any  attachment): 


fa)  Helium,  nitrogen  and  carbon  monoxide  ions  (and  electrons)  are 
generated  as  indicated  above. 

(b)  Helium  and  nitrogen  ions  charge  transfer  to  carbon  monoxide 
ions  infinitely  fast. 

(c)  Because  of  (b),  the  effective  generation  rate  for  the  carbon  monoxide 
ions  is  the  sum  of  all  the  generation  rates. 

+  +  _Q 

(d)  The  time  to  form  (CO)_  or  (CO)  is  estimated  to  be  about' 10  ‘  sec. 

«  X 

Since  this  time  is  short  compared  to  the  recombination  time,  we 
may  assume  that  this  process  proceeds  arbitrarily  fast. 

(e)  The  recombination,  of  (CO)*- or  (CO)  indicated  above  is  then  the 

t  x 

only  important  loos  process  for  electrons. 

Although  the  electron -ion  recombination  process  for  carbon  monoxide  (and 
other  ions)  involves  a  larger  molecule,  i.  e. ,  (CO)*  we  shall  simplify  our 


V£3J&Trv^>  t^'-'  e 
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notation,  whenever  possible,  by  writing  CO  instead  of  (CO/^  .  We 
shall  use  this  convention  for  other  ions  as  well. 

The  recombination  rate  for  carbon  monoxide  ions  is  i(C0+)  n  (CO  )  n  where 
+  ^  e 

a  (CO  )  is  the  recombination  rate  constant  for  CO  .  For  the  :ase  of  no 

+  +  2 
oxygen,  n(CO  )  =  ng,  and  the  recombination  rate  reduces  to  <r(CO  )ng  .  Using 

this  loss  teriivEqs.  (la)  and  (lb)  reduce  to 

dn  /dt  =  G(CO+)  +  GfHe1)  +  G(N*)  -  a(COt')n2  .  (2) 

e  6  6 

It  is  also  assumed  that  the  relative  concentration  of  oxygen  atoms,  when 

-2 

present,  is  less  than  about  10  ,  and  therefore  the  number  of  electrons 

derived  from  positive  oxygen  ionization  is  negligible.  The  effect  however, 
of  even  a  small  amount  of  oxygen  on  electron  loss  through  attachment,  may 
be  severe,  and  the  coupled  set  of  Eqs.  (la,  b  and  c)  must  then  be  solved. 

The  statements  (a)  and  (d)  are  uneffected  by  oxygen  and  they  still  serve  to 

simplify  Eq.  (lb).  In  Eq.  (Ic),  G(0~)  represents  the  attachment  process 

and  L(CL )  the  ion-ion  recombination  process.  We  shall  denote  the  oxygen 
c* 

recorribination  coefficient  with  a( O 2).  The  simplified  set  of  equations 
for  a  helium,  nitrogen,  carbon  monoxide  and  oxygen  mix  is: 


n  -  n(CO+)  -  n(0~)  (3a) 

e  c 

dn(CO+)/dt  ■:  G(CO+)  +  G(He+)  +  G(N*) 

-  <*(CO+)n  n(CO+)  -  a(Or)n(Or)n(CO+)  (3b) 

e  l  L 

dn(0")/dt  =  ng  0'  -  a(0“)n(0~)n(C0+)  (3c) 

where  (discussed  in  Section  3.  4)  is  the  attachment  rate  constant.  We 
have  ceased  distinguishing  between  O  and  (or  other  negative  oxygen 
ions). 

Because  of  the  very  anisotropic  electron  velocity  distribution,  it  is  convenient 

to  divide  the  electrons  into  two  classes,  as  indicated  earlier,  and  to  divide  the 

+  + 

ionization  function  G(j  )  into  two  parts.  Thus,  G(j  )  is  written  as 

G(j+)  =  in(j)g(j)  +  nen(j)  F(Te,  j)  (4) 

where  the  first  term  on  the  right  hand  side  corresponds  to  the  ionization  caused  by 
the  electron  beam  and  the  second  term  the  ionization  caused  by  the  electrons 
having  random  energies  of  the  order  of  1  eV.  n(j)  is  the  molecular  or  atomic 
number  density  of  species  j,  i  is  the  electron  beam  current,  and  g(j)  is  a 
normalized  ionization  rate  for  species  j.  F(Te»  j)  is  an  ionization  rate 
per  molecule  or  atom,  which  depends  on  the  electron  energy.  It  shall  be 
assumed  that  the  random  motion  of  the  electrons  can  be  approximated  by  a 
Boltzmann -Maxwell  distribution  at  temperature  Te.  The  functions  g(j) 
and  F(Te,  j)  are  discussed  in  Sections  3.1  and  3.2  respectively. 

g  g 

*  Although  Nighan  ’  has  calculated  that  the  electron  velocity  distribution 
in  CO2  and  CO  gas  mixtures  for  the  usual  discharge  conditions  will  depart 
from  a  Maxwellian,  we  will  often  treat  this  approximately  i ootropic  random 
motion  as  thermal. 
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2. 2  The  Analytic  Solution  of  Electron  Penalty  for  a  Single  Component 
Gas,  In  this  section  the  basic  equation  for  the  electron  density  (Eq.  1) 
will  be  simplified  by  considering  only  one  gas,  for  which  an  analytic  solution 
can  be  obtained.  An  approximate  equation  that  includes  the  effects  of  oxygen 
will  also  be  derived  and  solved. 

For  this  example  we  shall  assume  that  the  fill  gas  is  argon  with  a  small 
amount  of  oxygen.  Since  the  percentage  of  oxygen  is  assumed  small,  we 
neglect  the  positive  ions  of  oxygen.  Eqs.  (la,  b  and  c)  now  reduce  to: 


nQ  a  n(Ar  )  -  n(02) 


dn(Ar  )/dt  =  n(Ar)ig(Ar)  +  n  n(Ar)  F(Te,  Ar) 
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-  a(Ar+)n(Ar+)ne  -  a(C>2)n(02)n(Ar+) 
dn(0‘)/di  =  0‘  ne  -  a(0')n(02)n(Ar+)  . 


This  set  of  equations  cannot  be  solved  simply  without  numerical  methods. 

We  shall  obtain  analytical  solutions  for  three  cases:  (a)  no  thermal 
generation  (F(T)  *  0)  and  no  oxygen  attachment  (  0*  =  0),  (b)  thermal 
generation  but  no  attachment,  and  (c)  no  thermal  generation  but  attach¬ 
ment.  For  cases  (a)  and  (b),  we  obtain  simple  solutions  to  the  exact 
equations.  For  (c)  we  must  approximate  Eq.  (5)  in  order  to  solve  for  n#. 

We  shall  assume  that  the  electron  beam  is  on  for  a  length  of  time  p.  Figure  2 
shows  the  value  of  the  electron  density  as  a  function  of  time  for  the  cases  (a), 
(b)  and  (c). 

Case  (a):  For  F(l^)«  0  and  0'  *  0  we  have 

2 

dn  /dt  ■  -  an  +  c  (6a) 

o  e 

where  a  «  <*(Ar+)  and  c  *  n(Ar)ig(Ar) 
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Figure  2.  Electron  beam  current  and  electron  density  vs  time. 


The  solution  is: 

1  tO  •— 

For  the  buildup  of  n  :  n  =  (c/a)  '  tanh  Vac  t 

6  6 


o<  t<p 


For  the  decay  of  n^:  ng  =  (a(t  -  p)  +  1/n^)* 


where  n  =  n  (t  =  p'. 
eo  e 


Eq.  (6b)  and  (6c)  are  plotted  in  Figure  2. 


If  we  define  the  characteristic  rise  and  decay  times  to  be  T  *  (ac)” 

-1  r 
and  Tj  =  (n^a)  respectively,  then  we  note  that  =  Tj  if  the  build  up 

transient  is  complete  (i.  e. ,  if  n  ar  (c/a)^2).  A  sample  calculation 

for  760  torr  argon  with  1  mA/cm  electron  beam  current  yields  a  steady 

12  3 

state  equilibrium  electron  density,  n  ,  of  4.  3  x  10  /cm  and  a  rise  time 

®*l  m<j  j  jg 

of  1  H s.  This  was  calculated  using  a  =  8.  8  x  10~  cm  / a  and  c  =  1.2  x  10 
(cm  s)  . 


Case  (b]l:  For  F(T^  0  and  /?'  =  0  we  have 


dn  /dt  =  -a  n  -bn  +c 
e  e  e 


where  we  define  b  *  fi*  -  n(Ar)F(T  ),  In  this  case  b  simplifies  lo 

e 

«n(Ar)F(Te)»  The  solution  for  this  case  is: 


For  the  buildup  of  n  ;  n  *  —  -  ■ — 

("»b  +  R)  e  +  b  +  R 


Oit^p  (7b) 


For  the  decay  of  n  ;  ns 
0  © 


(!^»  i)eb<‘-p’ . 


where  R  «  (b2  +  4  ac)ly^2. 
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Case  (c):  For  this  case  Eqs.  (2a,  b,  and  c)  can  be  combined  to  yield 

dn  /dt  =-an2  (1  +  n(0")/n  )  -  bn  +  c.  (8) 

e  e  L  e  e 

Considering  attachment  and  nfChj/n  <c  1,  Eq.  (8)  is  of  the  form  of  (7a) 

c.  e 

and  its  solution  is  again  given  by  (7b)  and  (7c),  but  the  coefficient 
b  =  0'  for  this  case. 

The  value  of  n(0~)  'n  can  be  easily  estimated  for  steady  state  conditions, 
z  e 

From  Eq.  (5c)  one  calculates  - 

n(0~)  =  [0'/  a (0“)]  a  +n(0*)/ner1  (9) 

which,  for  n(0_)/n  small  compared  to  one,  gives  n(0~)  0'  /  a{0  ). 

“  e  3  1  .  —  ^ 

For  typical  values  of  0'  -  5  x  10  s  and  0(0,)  =  10”  cm^s"*,  one  calculates 

-9-3  12  -3 

n(C^)  =  5  x  10  cm  .  For  ng  =  10  cm  the  above  approximation  to 

Eq.  (8)  for  steady  state  conditions  is  well  justified.  During  the  initial 
build  up  of  n^,  the  ionization  term  c  is  more  important  than  the 
electron  recombination  term  and  therefore  the  approximation  to  Eq.  (8) 
should  also  be  good  for  this  period.  Late  in  the  after-glow^  the  recombina¬ 
tion  term  is  less  important  than  the  attachment  term  and  again  the  approxi¬ 
mation  to  Eq.  (8)  should  cause  no  significant  error  in  calculating  n  for 
this  period. 

If  the  turn  on  transient  is  complete  and  steady  state  is  achieved,  the 
equilibrium  electron  density,  n  for  cases  (b)  and  (c)  is  - 

n  =  [-b  +  (b2  +  4ac)1/2}  (2a)"1.  (10) 

eq 
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The  solution  for  n  for  cases  (a),  (b),  and  (c)  is  shown  in  Figure  2.  The 
e 

initial  rise  in  electron  density  is  governed  solely  by  the  electron  beam 
ionization  rate.  The  final  phase  of  the  electron  density  build-up  is  only 
slightly  influenced  by  attachment  and  thermal  ionization;  it  proceeds  faster 
for  the  case  of  attachment  and  slower  for  the  case  of  thermal  generation.  - 
The  rate  of  decay  of  the  electron  density  is  strongly  influenced  by  attachment 
and  thermal  generation  as  shown  in  Figure  2. 

2. 3  The  Electron  Density  Function  for  a  Multicomponent  Gas.  The 
simplified  basic  equations  (3a,  b  and  c)  for  a  gas  mixture  are  of  the  same 
form  as  those  for  a  single  component  gas  (5a,  b  and  c).  Therefore, 
the  solutions  for  n^  in  equations  (7b,  c)  and  (9),  also  represent  analytic 
solutions  for  the  multicomponent  gas  case,  provided  that  the  definition  of 
the  constants  a,  b,  and  c  is  generalized.  It  should  be  noted  that 
the  differential  equation  (7a)  is  approximate  when  electron  attachment 
is  significant,  since  n(j+)  is  not  equal  to  n^. 

For  the  multicomponent  gas,  the  coefficients  a,  b»  and  c  to  be  used  in 
the  approximate  solution  (7b,  c)  are: 

a  =  o(j+)  (11a) 

where  (j+)  is  the  final  ion  type  in  the  rapid  charge  transfer  process, 
e.g. ,  CO+  or  CO*  , 


b=  -  X  n(j)F(T,j)  (lib) 

j 

where  /?', derived  in  Section  3.  4,includes  the  2  body  and  3  body  attach¬ 
ment  process  and  n(j)F(Te,  j)  is  the  thermal  ionization  rate  of  atom 
species  j,  and 

c  =  i£n(j)g(j)  (11c) 

3 


where  in(j)g(j)  is  the  rate  of  ionization  produced  by  the  electron  beam  in 
species  j. 


The  qualitative  behavior*  of  n  is  depicted  in  Figure  3  for  various 
cases.  An  electron  beam  pulse  of  100  /is  duration  is  assumed  and 
results  for  pure  and  multicomponent  gases  with  and  without 
attachment  and  thermal  generation  are  shown.  In  the  next  section, 
numerical  solutions  of  Eqs.  (3a,  b,  c)  are  presented  and  discussed. 

In  the  section  following,  comparisons  are  made  between  the  calculated 
and  measured  results. 

2.4  Solutions  for  Various  Operating  Conditions.  Computer  solutions  of 
ne  for  various  operating  conditions  are  presented  in  tide  section.  A 
100  keV  electron  beam  pulse  of  100  ft s  duration  and  a  total  gas  pressure  of 
760  torr  at  273°K  was  assumed  for  each  case.  Table  I  lists  die  cases 
considered  with  parameters  chosen  to  illustrate  the  effect  on  the  electron 
density  caused  by: 

Changing  the  beam  current;  Figv  c 

Including  attachment;  Figures  4d,  6b  and  7b, 

Including  thermal  generation;  Figures  4e,  6c  and  7c, 

Changing  gas;  figures  5a,  b,  c,  and  d,  and 
Using  differer;  laser  gas  mixtures;  Figures  6  and  7. 

With  no  attachment  or  thermal  generation  the  effect  of  changing  the  beam 
current  is  as  anticipated:  the  maximum  electron  density  varies  as  the 
square  root  of  the  current,  and  the  rise  and  fall  times  decrease  at 
higher  currents.  From  Bq.  (4)  one  notes  that  a  change  in  the  gas 
density  is  equivalent  to  a  change  in  the  beam  current. 


The  electrom  beam 


0 


100  III 


For  pure  CO  with  no  thermal 
generation.  (High  generation 
and  recombination  rates) 


For  pure  He  with  no  thermal 
generation.  (Low  generation 
and  very  low  recombination 
rates.  ) 


For  CO,  N£,  and  He  with  no 
thermal  generation  and 
attachment. 


For  CO,  N2  and  He  with  thermal 
generation  but  no  attachment. 


For  CO,  N2,  and  He  with 
attachment  but  no  thermal 
generation. 


Figure  3.  Lffects  of  various  basic 


processes  on  electron  density. 


TABLE  1  -  Table  of  Rates  used  in  Figurea 
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Figure  4.  Electron  Density  for  Nitrogen  under  Various  Conditions 
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Ls  Figure  5  the  effect  of  the  low  generation  rate  and  the  very  low 

recombination  ratein  helium  is  noted;  namely  long  rise  and  fall  times 

12  3 

but  a  high  steady  state  electron  density  (its  value  would  be  7  x  10  /cm 
for  a  sufficiently  long  pulse).  The  formation  of  negative  ions  is  to  be 
expected  for  CO  and  CO2  discharges,  but  it  is  not  considered  in  Figures 
5c,  and  d.  - 


it' 


Figures  6  and  ?  give  the  electron  density  for  typical  laser  gas  mixtures. 


1 


Whenever  helium  is  present,  there  is  rapid  charge  transfer  to  carbon  ■ 
monoxide  or  carbon  dioxide.  Since  the  generation  rateof  helium  is  lew 
compared  to  N^,  CO  or  CQ2,  there  is  only  a  small  contribution  to  the 
electron  density  from  the  helium  (compare'Figures  6a  and  6d). 


The  purity  of  CO  is  often  a  problem  because  carbonyls  are  often  present 
and  tney  attach  electrons  rapidly.  ^  The  inclusion  of  some  oxygen  attach¬ 
ment  should  approximate  the  effect  of  carbonyl  electron  attachment.  Little 
is  known  about  the  recombination  of  carbon  dioxide  because,  in  most 
discharges,  many  of  the  free  electrons  attach  and  form  negative  ions  which 
mask  the  recombination  process.  (CO^  has  a  large  cross  section  for  negative 
ion  formation  at  electron  energies  of  about  4  eV). 


fc 

£ 


I 

t 


2.5  Comparison  of  Measured  and  Calculated  Currents.  In  electrical 
discharges,  operating  at  atmospheric  pressure  with  a  low  degree  of  ioniza¬ 
tion,  it  is  difficult  to  measure  the  electron  density  directly  (by  probes 
or  by  microwave  methods).  A  measure  of  the  accuracy  of  the  electron 
density  calculations  may  be  obtained  by  comparing  the  current  which  flows 
between  the  sustainer  electrodes  in  the  laser  plasma  region  (see  Vigure  1). 
The  expected  value  of  current  density  i  is 

*  I  vd  I  <U’ 
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where  e  is  the  electron  charge  and  v^  the  electron  drift  velocity 

in  the  electric  field  of  the  sustainer  electrodes.  The  values  for  vj 

12 

were  obtained  from  Brown,  and  the  electric  field  was  taken  to  be 
the  voltage  between  the  electrodes  divided  by  the  separation.  The 
total  measured  sustainer  current  divided  by  the  area  gave  the  average 
current  density  i&. 

The  results  of  preliminary  measurements  of  discharge  currents  in  nitrogen 

and  in  carbon  monoxide  are  shown  in  Figures  8  and  9.  The  electron  beam 

2 

current  density  was  about  0.  5  mA/cm  ,  and  the  current  measurements  in 
the  laser  cavity  were  made  after  the  electron  density  build  up  was 
complete.  Both  ft'  and  F(Tg)  were  taken  to  be  zero  in  the  calculation. 

The  cathode  fall,  estimated  to  be  about  400  volts,  was  neglected. 

The  agreement  between  the  measured  and  calculated  discharge  current 
values  for  nitrogen  was  better  than  expected  with  this  crude  model.  The 
lower  measured  current  values  for  carbon  monoxide  may  be  due  to  some 
attachment  process. 


Measured 
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3.  0  c  BASIC  PHYSICAL  PROCESSES 


In  this  section  we  shall  discuss  the  basic  reactions  and  the  reaction 
rates  giving  rise  to: 

(1)  Ionization  of  the  fill  gas  atoms  by  the  electrons  of  the  E-beam, 

(2)  Ionization  of  the  atoms  by  electrons  having  about  1  eV  thermal 
energy, 

Electron -ion  and  ion -ion  recombination. 

Electron  attachment  to  and  detachment  from  oxygen,  and 
Charge  transfer  processes. 


(3) 

(4) 

(5) 


3. 1  Ionization  of  Fill  Gas  Atoms  by  Electrons  of  the  E-Beam.  The 


reactions, by  which  the  very  energetic  electrons  from  the  electron  gun 
ionize  neutral  gas  particles,  are  cascade  processes,  which  may  be 
expressed  as  follows  using  Ar. 


e(E)  +  Ar  -*Ar  +  e(E  -  V.  -  Wj)  +  e(W1) 


e(W1)  +  Ar^-Ar  +  e(W1  -  V  -  Wy  +  e(W2) 


(13a) 

(13b) 


where  E  is  the  energy  of  the  beam  electron  (E-100  keV),  V.  is  the  ioniza* 


tion  energy  and  W  is  the  kinetic  energy  of  the  ejected  electron. 


The  energetic  electron  ionizes  an  argon  atom  and  looses  an  mount  of 
energy  approximately  equal  to  the  ionization  energy  plus  the  kinetic  energy, 

Wj,  of  the  secondary  electron.  The  cross  sections  <r.  for  this  process, 
listed  in  Table  II,  are  discussed  below.  The  generation  process,  however, 
is  not  complete  as  the  secondary  electron  at  energy  may  then  ionize 
additional  atoms.  The  likelihood  of  this  latter  process  is  more  difficult 
to  calculate.  The  highly  energetic  electron  may  also  impart  some  further 
energy  to  the  argon  atom,  however,  this  is  not  included  in  the  above  equations. 


u 
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Ye  will  first  discuss  the  cross  section  catenation  for  process  (13a)  and 
the  distribrikB  of  kinetic  energy  (Wj)  of  tte  ionised  electron,  and  then 
compare  these  energy  losses  of  the  primary  electa  on  with  stopping  power 
data.  We  will  then  discuss  process  J13b)  and  denote  with 
<  the  ""mh»r  q£  forfoer  ionisations  caused  by  the.  electron  at  energy  Vr 
Finally,  we  will  relate  the  total  number  of  ion  pairs/cm  caused  by  the  passage 
of  one  high  energy  electron  to  the  results  obtained  by  a  rule  of  thumb 
calculation.  The  latter  is  based  on  calculating  the  ion  pair  generation  rate 
from  the  stopping  power  by  assuming  that  each  ion  pair  takes  about  30  -  35  eV 
energy  of  the  primary  electron. 


51 

I 


The  cross  sections  for  process  (13a)  base  been  measured  by  many  investigators 
at  low  energies  (  <  1  keV).  Less  work  has  been  done  at  energies  around 
100  keV.  Sc  hr  am  et  al**  have  made  measurements  of  foe  first  eight  gases 
listed  in  Table  II  at  energies  up  to  20  keV  and  have  tabulated  the  constants 
necessary  to  nse  the  Bethe  formula}^  which  allows  one  to  calculate  the  ionisa¬ 
tion  cross  sections  for  100  keV  primary  electrons  (or  other  energies  of 
interest). 


Using  foe  very  simple  colli sional  model  presented  in  reference  15,  (page 
276),  winch  considers  a  highly  energetic  electron  encountering  other 
electrons  at  rest  (process  13a),  one  can  calculate  ff.  and  foe  average 
kinetic  energy  of  foe  ejected  electron  <  Wj  >  .  For  E  =  100  keV  and 
Vj  =  15  eV,  one  calculates  =  0. 043  x  10’18Z  cm2  and  <  Wj  >  *  105  eV, 
where  Z  is  foe  atomic  (or  molecular)  number.  The  calculated  value  of 
S-  is  about  a  factor  of  turn  too  small  when  compared  with  measured  values  in 
Table  II.  The  value  of  <  Wj  >  is  within  about  20%  of  foe  value  one  obtains 
from  stopping  power  data.  For  example,  in  argon  at  1  atmosphere,  with 
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9.  =  L  9  x  10  cm  .  there  are  51  direct  iororiag  processes  per  cm 


for  each  primary  electron.  The  energy  loss  per  cm  for  each  primary 


electron  is  thus  51  (Vj  4  <  Wj  >  )  =  6. 1  ktV/cm.  The  tabulated  energy 
loss  given  by  Berger  and  Seltxer*^  is  5. 1  kV/cm. 


The  kinetic  energy  of  the  ejected  electron  decays  from  a  value  of  W 
toward  a  steady  state  value  (—  1  eVj.  There  are  many  processes  which 
compete  for  this  excess  kinetic  energy  and  we  will  not  attempt  to  discuss 


how  Bus  energy  is  split  between  ionising  and  nonionising  processes. 
17 


Fano  proposed  a  theory  which  states  that  the  average  energy*  w,  which 
is  expended  in  prodiicing  an  ion  pair  by  foe  processes  of  Eq.  (13a)  and  (13b), 
and  subsequent  ionization  processes,  should  be  quite  independent  of  foe  atom 
structure  and  is  about  30  eV.  This  is  born  out  experimentally  and  more 
precise  values,  taken  from  reference  15  (page  233).  are  given  below.  Also 
shown  below  are  some  measured  values  of  £ ,  foe  number  of  further  ion 
pairs  produced  by  foe  electron  ejected  in  foe  primary  ionization  process, 

Bq.  (13a). 

w(ev) 
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Hydrogen 

Nitrogen 

Argon 

Oxygen 

Helium 

Neon 

Krypton 

Xenon 

CO 


CO, 


35 

36.2 
26 

32.3 
30 
28 
24 
22 
34 
34 


1.3 

2.3 
2.7 
2.5 


The  data  used  in  the  generation  rate  calculations  is  presented  in  Table  II. 

2 

The  ion  pair  generation  rate,  g,  is  normalized  for  1  mA/cm  primary 

3 

electron  beam  current  at  100  kev  and  for  1  gas  atom  or  molecule/cm  . 


4 

3 


t<5 


It  «u  calculated  by  dividing  the  value  of  the  stopping  power  of  the  gas 
taken  few  Berger*^  by  w,  the  average  energy  loss  per  ion  pair.  Also 
listed  are  the  stopping  powers*”  and  the  ionisation  cross  sections  for 
the  production  of  ion  pairs  directly  by  the  primary  electrons.  The 
generation  rates  in  Table  I  were  obtained  by  multiplying  g  by  the  beam 
current  (in  mA)  and  the  fill  gas  number  density. 


3. 2.  Ionisation  by  Electrons  having  about  1  eV  Thermal  Eneri 


electrons  in  .the  plasma,  produced  by  the  ionising  processes  described 
above,  are  maintained  at  an  average  energy  of  about  1  eV.  The  motion  : 
of  these  electrons  is  primarily  random  with  a  relatively  slow  drift 
( —10  cm/s)  toward  the  positive  electrode.  Because  of  the  pre¬ 
ponderance  of  electron -neutral  collisions  (compared  to  electron -electron 
collisions)  and  the  non -thermal  conditions  (gas  kinetic  energy  —  0.  03  eV), 
the  electron  velocity  distribution  is  non -Maxwellian  (see  References  8  and  9). 


The  correct  method  for  the  calculation  of  the  ionization  rate  caused  by  these 
plasma  electrons  would  be  to  use  the  velocity  distribution  appropriate  to  the 
particular  gas  mixture  and  the  applied  electric  field.  A  much  simpler 
procedure  used  here,  which  should  yield  rates  within  a  factor  of  about 
two,  is  to  take  the  Maxwellian  velocity  distribution  equivalent  in  temperature 
to  2/3  of  the  average  electron  energv.  The  equation  for  F(Te,  j)  is 

co 

F(Te,  j)  =  J  a(j,  s)  s  f(s)  ds  (14) 

o 

where  j  denotes  the  type  of  atom  or  molecule,  s  and  f(s)  are  the 

electron  speed  and  speed  distribution  function  and  a  (j,  s)  is  the  ionization 

cross  section.  The  values  for  a  were  taken  from  D.  Rapp  and  P.  Englander - 
19 

Golden.  The  results  of  this  integration  are  plotted  in  Figure  10  as  a  function 
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of  electron  temperature.  In  the  above  treatment  only  direct  two 
body  ionizing  collisions  are  considered  and  ionization  via  meta  stable 
states  or  other  processes  are  ignored. 


3. 3  .Electron -Ion  and  Ion -Ion  Recombination.  The  electron-ion  recombina¬ 

tion  process  represents  fee  reverse  of  the  electron  impact  ionization  process. 
The  ion-ion  recombination  process  is  responsible  for  the  neutralization  of  die 
negative  ions  generated  by  the  electron  attachment  process.  The  recombination 
processes  are  characterized  by  a  rate  coefficient  a.  The  number  of  electron- 

o  •  « 

ion  recombination  events  per  crri  sec  is  given  by  a(j  )n(j  jn^,  where  a(j+) 
depends  on  the  ion  species.  For  ion -ion  recombination,  O  depends  on  both  ion 
species.  For  example,  for  <>£  and  ,  the  recombination  rate  is 
«(02.  N2)n(02)n(N2).  These  two  different  recombination  processes  will  be 
discussed  separately. 


Of  the  many  different  types  of  electron-ion  recombination  processes 

20 

possible  .3ee,  for  example,  McDaniel  ),  it  appears  that,  at  gas  pressures 
of  the  order  of  1  atmosphere,  dissociative  recombination  dominates.  In 
particular,  for  diatomic  molecules  for  example, 


+  N,  N*  ,  then  N*  +  e  2N_ 
C  C  4  4  C 


(15) 


and 

CO+ +  n(CO)^(CO)+11,  then  (CO)+.,  +  e-*(n  +  1)  CO  (16) 

Ti+i  n+x 

seem  to  be  the  important  processes.^’  ^  For  the  noble  gases, a  similar 
process  accounts  for  the  recombination, 


X+  +  X  -#-X*  ,  then  X*  +  e  -*.2X, 


(17) 
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'where  X  represents  one  of  the  noble  gas  atoms.  The  dependence  of  some  of  the 

10 

recombination  processes  on  electron  energy  have  been  observed  and  the 
values  quoted  in  Table- III  are  for  the  temperatures  indicated.  The 
author  is  not  aware  of  a  value  of  a  for  electron-carbon  dioxide 
recombination.  In  order  to  make  the  calculations^  we  have  used  the 
measured  value  of  a  for  electron-carbon  monoxide  recombination. 

References  to  numerous  other  measurements  of  a  are  given  by  Oskam 

24  iZ 

and  Mittelstadt  and  Hasted. 


The  author  is  not  aware  of  ion -ion  recombination  coefficients  for  negative 

22 

oxygen  ions  and  positive  carbon  monoxide  or  carbon  dioxide  ions.  Sayers 
has  measured  a  for  oxygen  and  nitrogen.  Because  of  the  similarity  between 
CO  and  ^  in  the  ionization  and  electron  recombination  processes  we  have 
assumed  an  equal  recombination  coefficient.  For  the  lack  of  a  better 
value  we  have  also  assumed  the  same  a  for  CO^. 

3.4  Electron  Attachment  to  and  Detachment  from  Oxygen.  The  electrons 
in  the  discharge  have  a  propensity  to  attach  themselves  to  any  oxygen  molecules 
present  in  the  system.  In  the  process  of  attachment,  the  electron  and 
oxygen  molecule  evolve  energy  (approximately  1.  5  eY  for  O  and  0.  4  eV  for 
O”)  which  has  to  be  removed  some  way.  We  would  then  expect  to  have  a 

M 

variety  of  processes  for  attachment  (as  is  the  case  for  recombination). 
Phenomenologically  these  processes  are  divided  into  two  classes:  two 
body  attachment  and  three  body  attachment.  One  may  then  write  the  rate 
of  electron  attachment  as: 


B'  n 
H  e 


nn(OJ[^  Z  K(/)n(i)] 
e  2  i 


(18) 


where  the  /  represents  a  particle  species  and  the  sum  is  to  be  extended 
over  all  the  species  present.  fi  and  K  are  the  two  and  three  body  attach¬ 
ment  coefficients  respectively.  These  coefficients  depend  strongly  on 
electron  temperature  as  discu38td  in  references  11  and  26.  For  a  gas 
mixture  of  helium,  nitrogen,  carbon  monoxide,  and  oxygen  Eq.  (18)  becomes 
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TABLE  ID 


Electron -loo  Recombination  Coefficicats 


Gas 

a  (cm3/*) 

Remarks 

«•* 

-9 

4  -  10  7 

At -50  torr  and  T  at.  63  eV24 

N'2 

3.4-  10~7 

At -20  torr  and  T  — .03  eV° 

.  ♦ 

AT 

-7 

8.8-10 

At  —20  torr  atd  T-.U  eV*3 

-4 

1.2*10~& 

A* —20  torr  aadTi.WtV24 

1. 4*  10*^ 

Af°20  torr  and  T  a=.  03  eV“4 

< 

3  '  10~8 

10 

At  l  atm  and  T  —I  eV 

<ccC 

1  -io~7 

10 

At  1  atm  and  T  —  1  e  V 

Ion-Ion  Recombination  Coefficient 


-  1*  10 


6 


Air 


At  1  atm 
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I1  *  V*V  *  EjHrNHrj 


T&e  «zlux  «f  ^  »wl  ta.  nforts  4L  ft  u£  7  ^»ge  oBoAhi  fintt  £%.  fcHI 


*«£3  J 

#  =  2  at  SB  can  Is,  crfr  apsfiaterf  ftm  FSpcro  $  of  ttefiecencc  7ft 

SerT  =IeT 
c 

KgO.^  sLSeM  ^cas.^  Is.  esfiraasSated  frooa  F%ut  E  off  BUfitratct  2ft 

-Ssy  Te  =  2  eV 

EfcKal  =  Eff^^cas*  /*_  extrapolated  frwx  Fiptre  E  off  Btefferetoce  2ft 

Ssrr  Te  »  £  eV 

_jj  ft 

KgX^D  =  %  *  B®  cm  f s,  eCrxjaTiitrc  €rooa  Ftgysc  12  of  Reference  2ft 

Gar  Tts!*V. 

TESte  "ft""  va&ces  Gar  CO  aai  QG>  are  cAapm  So  t&e  arfftor  asnff  sac  taerefiiore 

WKBXaX^  SCftC. 


-ss.  ULBacbr 


Tke  reverse 


agyar  to  gcoceec  at  a  very  Saw  rate  at  a 

13  5 

and  at  Crrrsfties  off  0,?  less  ftas.  19  Icre  .so  ftat  rt 
Tfce  oetac&ssesl  rate  at  0§*K  is  I*  ** 


°2“C2 


less  ftatr 


sfijWOjS. 


;e  411. 


Tie  attachment  of  electrons  to  CO  agpearo  to  fte  aarimpsrtaaat  at  ioxr  electron 
temper aftcrcs.  At  tie  preseat,  tee  utter's  oaty  £*  at 

Edercace!!.  eitdi  does  met  tree  fata  fee  electron  «erpe»  Idee  f  eV. 

Tfce  attwrWejit  rate  for  CO*.  Severn,  proceeds  frits  fast  at  electron  energies 
off  4  eV.  ** 
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3.5  Ctett  Ttawfeg  Processes.  Eteftenmc  dar|e  trmfer  processes  | 

aaaai  ftc  sobs  oi  mokocles  a  i  gas  £<dor|e  proceed  npifif.  The  rate 
contents,  k.  for  eozr  of  the  expected  processes  have  been  measured 
asi  are  fhta  la  Table  IV.  hr  order  to  help  sort  out  fire  rates  of  these  pro¬ 
cesses  «e  Medicate  a  characteristic  lifetime  for  the  ion  which  appears  on 
the  left  side  of  the  react  sort  tipitian.  For  example,  the  rate  equation  for 
fie  hrifozn-asfcrogent  reaction1  is 

«(BeV*  =  -ha(He^J  *{^1.  ~ 

Ve  define  fee  characteristic  lifetime  for  fire  helium  ion  to  be  (k  a(Nj)j'V 

The  characteristic  times  given  in  Table  £V  assume  a  total  gas  pressure 

0. 

of  763  flarr  at  273  K_  and  a  mixture  ratio  of  5»7/l  for  He.  and  CO 
for  CO>r  ana  £  burr  af  Oj. 

Chrster  fomalxB  rates  for  foe  Xl  molecule  are  also  listed  in  Table  IV, 
and  for  cwwparisaa.  some  recombination  rates  for  electron-ion  and  ion- 
foe  processes  are  tachadedL  For  foe  electron-ion  recombination 
characteristic  times,  an  electron  density  of  10^  /cm^  is  assumed  and, 
for  foe  tsa-iae  process  a  nitrogen  ion  density  of  10^  /cm*  is  assumed. 

The.  identity  of  foe  recombining  oxygen  and  nitrogen  ions  is  not  known 
isc  therefore  foe-  final  equation  of  Table  IV  represents  foe  process 
schematically  only. 


The  cbaracteaistic  lifetimes  listed  in  Table  IV  show  that  foe  main  ionic 
charge  transfer  processes  proceed  very  Cast  compared  to  the  ion -electron 
t  ecensbiaatfoa  processes. 
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4. 0  SUMMARY 


In  an  electrical  discharge  gas  laser,  the  two  most  important  plasma 
parameters  for  electron  pumping  are  electron  density  and  electron 
energy.  In  this  report  the  important  physical  processes  which  control 
the  electron  density  are  discussed  and  a  set  of  equations  for  calculating 
the  density  are  formulated  and  solved.  It  is  shown  that  the  equilibrium 
value  of  the  electron  density  of  an  electron  beam  stabilized  CO  or  CO^ 
laser  discharge  may  be  calculated  approximately  from  a  simple  quadratic 
equation.  Because  of  rapid  charge  transfer,  the  only  significant  number 
of  positive  ions  present  in  the  discharge  are  those  of  CO  or  CO£.  The 
electron  generation  rate  is  the  sum  of  the  rates  for  all  the  various  gas 
constituents.  These  generation  rates  are  proportional  to  the  electron 
beam  current  density,  and  mass  density  of  the  fill  gas.  The  calculated 
electron  densities  from  these  formulations  are  in  reasonable  agreement 
with  measured  values. 
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